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ABSTRACT 

We present an unsupervised machine learning approach that can be employed for estimating photomet- 
ric redshifts. The proposed method is based on a vector quantization approach called Self-Organizing 
Mapping (SOM). A variety of photometrically derived input values were utilized from the Sloan Dig- 
ital Sky Survey's Main Galaxy Sample, Luminous Red Galaxy, and Quasar samples along with the 
PHATO data set from the PHoto-z Accuracy Testing project. Regression results obtained with this 
new approach were evaluated in terms of root mean square error (RMSE) to estimate the accuracy 
of the photometric redshift estimates. The results demonstrate competitive RMSE and outlier per- 
centages when compared with several other popular approaches such as Artificial Neural Networks 
and Gaussian Process Regression. SOM RMSE-results (using Az=z p / lot -z spec ) for the Main Galaxy 
Sample are 0.023, for the Luminous Red Galaxy sample 0.027, Quasars are 0.418, and PHATO syn- 
thetic data are 0.022. The results demonstrate that there are non-unique solutions for estimating 
SOM RMSEs. Further research is needed in order to find more robust estimation techniques using 
SOMs, but the results herein are a positive indication of their capabilities when compared with other 
well-known methods. 

Subject headings: methods: data analysis, methods: statistical, galaxies: distances and redshifts 



1. INTRODUCTION 

There is a pressing need for accurate estimates of 
galaxy photometric redshifts (photo-z's) as demon- 
strated by the increasing number of papers on this topic 
and especially by recent attempts to objectively com- 
pare methods (e.g. Hildebrandt et al. 2010; Abdalla et 
al. 2011). The need for photo-z's will only increase as 
larger and deeper surveys such as Pan-STARRS 5 (Kaiser 

2004) , LSST 6 (Ivezic et al. 2008) and Euclid (Sorba k 
Sawicki 2011) come on-line in the coming decade. The 
photometric-only surveys (Pan-STARRS, LSST) will 
have relatively small numbers of follow-up spectroscopic 
redshifts and will rely upon either template-fitting meth- 
ods such as Bayesian Photo-z's (Benftez 2000) Le Phare 
(Ilbert et al. 2006), or training-set methods such as those 
discussed herein. The Euclid mission may include a slit- 
less spectrograph offering far more training-set galaxies. 

A diverse set of regression techniques using training- 
set methods have been applied to the problem of estimat- 
ing photometric redshifts in the past 10 years. These 
include Artificial Neural Networks (Firth et al. 2003; 
Tagliaferri et al. 2003; Ball et al. 2004; Collister & Lahav 
2004; Vanzella et al. 2004), Decision Trees (Suchkov et al. 

2005) , Gaussian Process Regression (Way & Srivastava 
2006; Foster et al. 2009; Way et al. 2009; Bonfield et al. 
2010; Way 2011), Support Vector Machines (Wadadekar 
2005), Ensemble Modeling (Way et al. 2009), Random 
Forests Carliles et al. (2008), and Kd-Trees (Csabai et 
al. 2003) to name but a few. 
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On the other hand, even though Self-Organizing Maps 
(SOMs) have been used extensively in a number of other 
scientific fields (the paper that opened the field, Koho- 
nen (1982), currently has over 2000 citations) they have 
been used sparingly thus far in Astronomy (e.g. Mahdi 
2011; Nairn et al. 1997; Way, Gazis & Scargle 2011), and 
only this year in estimating photometric redshifts (Geach 
2011). 

In this work we attempt to use SOMs to estimate pho- 
tometric redshifts for several Sloan Digital Sky Survey 
(SDSS, York et al. 2000) derived catalogs of different 
galaxy types, including Quasars along with the PHATO 
data set of Hildebrandt et al. (2010). In Section 2 we 
describe the input data sets used, in Section 3 we give 
an overview of SOMs, and some conclusions in Section 
4. 

2. DATA 

Three different data sets derived from the SDSS Data 
Release Seven (DR7, Abazajian et al. 2009) were used. 
They include the Main Galaxy Sample (MGS, Strauss et 
al. 2002) the Luminous Red Galaxy Sample (LRG, Eisen- 
stein et al. 2001), and the Quasar sample (QSO, Schnei- 
der et al. 2007). Data from the Galaxy Zoo 7 (Lintott 
et al. 2008) Data Release 1 (Lintott et al. 2011) survey 
results were used to segregate galaxies as Spiral or Ellip- 
tical in the case of the MGS and LRG samples. Details of 
how this was done are given in Way (2011). Dereddened 
magnitudes (u,g,r,i,z) were used as inputs in all scenar- 
ios. The same SDSS photometric and redshift quality 
flags on the input variables were used as in Way (2011). 
In addition we used the simulation-based PHATO data 
set (see Hildebrandt et al. 2010) which was constructed 
to to test a variety of different photo-z estimation meth- 
ods. The PHATO data set consists of 5 SDSS like filters 
(u,g,r,i,z) used on MEGACAM at CFHT (Boulade et al. 

7 http://www.galaxyzoo.org 
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2003) with an additional 6 input niters (Y,J,H,K,Spitzer 
IRAC [3.6], Spitzer IRAC [4.5]) giving a total of 11 filters 
spanning a range of 4000A to 50,000A. This large range 
should help to avoid color-redshift degeneracies that can 
occur if ultraviolet or infrared bandpasses are not used 
(Benftez 2000). The PHAT0 synthetic photometry was 
created from the Le Phare photo-z code (Arnouts et al. 
2002; Ilbcrt ct al. 2006). Initially Le Phare creates noise 
free data, but given the desire to test more real-world 
conditions we utilized the PHAT0 data with added noise. 
A parametric form was used for the signal-to-noise as a 
function of magnitude where it acts as an exponential at 
fainter magnitudes and a power-law a brighter ones. The 
magnitude cut between these two regimes is filter depen- 
dent and is given in Table 2 of Hildcbrandt et al. (2010). 
The larger of two catalogs was used herein (as suggested 
for training-set methods) that contains ~ 170,000 ob- 
jects. 

Since we use a training-set method our original data 
sets are split into training=89%, testing=10% and val- 
idation=l%. Validation was only used in the Artificial 
Neural Network algorithm discussed in the next section. 
The full size of each input data set are listed in paren- 
theses in column 1 of Table 1. 

3. METHODS 

Several methods in use for calculating photometric red- 
shifts were compared with the SOM results: the Artifi- 
cial Neural Network code of Collister & Lahav (2004) 
(ANNz), the Gaussian Process Regression code of Fos- 
ter et al. (2009) (GPR), as well as simple Linear and 
Quadratic regression. The latter is comparable to that 
of the Polynomial fits used by Li & Yee (2008). Both the 
ANNz and GPR codes are freely downloadable 8 . Details 
on the ANNz and GPR algorithms can be found in their 
respective citations above. 

The main purpose of Self-Organized mapping is the 
ability of SOMs to transform a feature vector of ar- 
bitrary dimension drawn from the given feature space 
of photometric inputs (e.g., the SDSS u,g,r,i,z magni- 
tudes) into simplified 1- or 2-dimensional discrete maps. 
The method was originally developed by Kohonen (1982, 
2001) to organize information in a logical manner. This 
type of machine learning utilizes an unsupervised learn- 
ing scheme of vector quantization, known as competi- 
tive learning in the field of neural information process- 
ing. It is useful for analyzing complex data with a-priori 
unknown relationships that are visualized by the self- 
organization process (Kohonen 2001). 

A SOM is structured in two layers: an input layer and 
a Kohonen layer (Figure 1). For example, the Koho- 
nen layer could represent a structure with a single 2- 
dimensional map (lattice) consisting of neurons arranged 
in rows and columns. Each neuron of this discrete lat- 
tice is fixed and is fully connected with all source neu- 
rons in the input layer. For the given task of estimating 
photometric redshifts, a 5-dimensional feature vector of 
the u,g,r,i,z magnitudes is defined. One feature vector 
(u,g,r,i,z) is presented to 5 input layer neurons. This 
typically activates (stimulates) one neuron in the Koho- 
nen layer. Learning occurs during the self-organizing 

8 GPR: http://dashlink.arc.nasa.gov/algorithm/stableGP and 
ANNz: http://www.star.ucl.ac.uk/ lahav/annz.html 



Table 1. Results 



Data a Method b vrmse c Outlier d 

50% 10% 90% 



MGS GPR 0.02087 0.02072 0.02096 0.11629 

(455803) ANNz 0.02044 0.14482 

SOM 0.02339 0.1689 

Linear 0.02742 0.02729 0.02758 0.35986 

- Quadratic 0.02494 0.02412 0.02762 0.29184 

LRG GPR 0.02278 0.02256 0.02309 0.41898 

(143221) ANNz 0.02138 0.41176 

SOM 0.02689 0.64292 

Linear 0.02896 0.02896 0.02897 0.71516 

Quadratic 0.02382 0.02376 0.02402 0.45510 

MGS-ELL GPTT 0.01455 0.01434 0.01473 0.06591 

(45521) ANNz 0.01442 0.06591 

SOM 0.02044 0.10984 

Linear 0.01745 0.01731 0.01756 0.19772 

- Quadratic 0.01612 0.01609 0.01629 0.10984 

MGS-SP GPR 0.02078 0.02061 0.02093 0.13305 

(120266) ANNz 0.01991 0.05821 

SOM 0.02426 0.04158 

Linear 0.02539 0.02529 0.02555 0.28272 

- Quadratic 0.02326 0.02296 0.02607 0.20788 

LRG-SP GPR 0.01416 0.01397 0.01436 0.00000 

(13708) ANNz 0.01516 0.00000 

SOM 0.01848 0.07299 

Linear 0.01635 0.01627 0.01649 0.07299 

Quadratic 0.01469 0.01462 0.01477 0.00000 

LRG-ELL (JFK 0.01186 0.01162 0.01224 0.00000 

(27378) ANNz 0.01298 0.10961 

SOM 0.01568 0.00000 

Linear 0.01362 0.01361 0.01364 0.10961 

- Quadratic 0.01263 0.01254 0.01274 0.07307 

TJSO GTTT 0.37342 0.03967 0.37626 50.96627 

(56923) ANNz 0.65802 88.54533 

SOM 0.41821 54.23401 

Linear 0.57061 0.57010 0.57102 84.64512 

- Quadratic 0.53972 0.53679 0.54539 81.27196 

phatO GPR 0.01487 0.01436 0.01532 0.03539 

(169520) ANN 0.01805 0.05309 

SOM 0.02236 0.37754 

Linear 0.08703 0.08702 0.08704 19.34875 

Quadratic 0.02436 0.02433 0.02438 0.19467 



a MGS=Main Galaxy Sample (Strauss et al. 2002), 
LRG=Luminous Red Galaxies (Eisenstcin ct al. 2001), 
SP=Classificd as spiral by Galaxy Zoo, ELL=Classified as 
elliptical by Galaxy Zoo, QSO=Quasar sample (Schneider ct al. 
2007) 

b GPR=Gaussian Process Regression (Foster ct al. 2009), 
ANNz=Artificial Neural Network (Collister & Lahav 2004), 
SOM=Self-Organizing Maps (SOM-4100 and SOM-5100 see Fig- 
ure 2 for details), phatO=PHAT synthetic sample 

c We quote the bootstrapped 50%, 10% and 90% confidence levels 
as in Way et al. (2009) for the root mean square error (RMSE) 
when available. 

d Percentage of points defined as outliers. Following a prescription 
similar to that of Hildebrandt et al. (2010) we quote the percentage 
of points outside of Az=z p f lof — z spec ± 0.1 



procedure as feature vectors drawn from a training data 
set are presented to the input layer of the SOM network 
(Figure la). These feature vectors are also referred to as 
input vectors. Neurons of the Kohonen layer compete to 
see which neuron will be activated by the weight vectors 
that connect the input neurons and Kohonen neurons. In 
other words, the weight vectors identify which input vec- 
tor can represented by a single Kohonen neuron. Hence, 
they are used to determine only one activated neuron in 
the Kohonen layer after the winner-takes-all principle 
(Figure lb). 
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The SOM is considered as trained after learning, at 
which time the weights of the neurons have stored the 
inter-relations of all 5-dimensional u,g,r,i,z feature vec- 
tors. Then, known spectroscopic redshift values for all 
input vectors of a test data set that are represented by a 
single Kohonen neuron are averaged (Fig. lb). The red- 
shift mean value represents all 5-D u,g,r,i,z vectors that 
are similar to the weight vector of the activated Koho- 
nen neuron. The more Kohonen neurons there are the 
more precisely each input vector can be represented by 
a weight vector. However, the total number of Kohenen 
neurons are optimized for each data set (see Figure 2). A 
practical overview about the learning/training process is 
described by Klose (2006); Klose et al. (2008, 2010) and 
in much greater detail by Kohonen (2001). 

After training, the u,g,r,i,z input vectors of a test data 
set are presented to a trained SOM. At the end of a clas- 
sification step, every Kohonen neuron approximates an 
input vector whereby similar/dissimilar input data were 
represented by neighboring/distant neurons. One neuron 
could even classify several input vectors, if these input 
vectors were very similar compared to the other given 
input vectors. Results from the photometric redshift ap- 
proximations are then compared to known spectroscopic 
redshift data. Regression performance is estimated based 
on the root mean square error (RMSE) of the predicted 
photometric redshifts and the known spectroscopic red- 
shifts (using Az—Zphot-^spec)- To reiterate, during the 
training phase, each Kohonen neuron identifies a cer- 
tain number of galaxies that are characterized by similar 
u,g,r,i,z intensities. Photometric redshift data were then 
averaged for these intensity values. 

The SOM approximates the input feature space and 
maps it into an output space. The output space shows 
the SOM approximation as a 2-D map (Haykin 2009). 
Best results can be obtained with an optimization scheme 
such that the RMSE of the test data set is minimal as 
illustrated in Figure 2. Accuracy (e.g. RMSE) depends 
on the size of the Kohonen map. The number of neurons 
in the Kohonen map can be considered a regularization 
parameter (£) as shown in Figure 2. 

Figure 2 shows that RMSE is high when the number 
of Kohonen neurons is too small (£ <2000) or too large 
(£ > 10000) and hence that the 5-dimensional u,g,r,i,z- 
input space is underfit or ovcrfit. Theoretically, a global 
minimum of the RMSE-curve might exist. However, the 
input feature space for the given photometric redshift 
problem shows a very rough RMSE-curve (Figure 2) with 
at least 2 local minima. In this case it is clear that SDSS 
redshift estimation tends to have several local minima, 
which makes is important to chose the right optimiza- 
tion method to determine the SOM network size. On the 
other hand, the smoother the RMSE-curve is the better 
gradient methods can be utilized. Evolution strategies or 
genetic programming could be applied to rougher curves 
with many local minima. This in turn can make it cum- 
bersome to find fast back-propagation Artificial Neural 
Network (ANN) structures, especially when data sets are 
small. 

Another advantage of SOMs in comparison to ANNs is 
that there is no need to optimize the structure of SOMs 
(e.g., number of hidden layers), since it is based on un- 
supervised learning. 

Only the size of the Kohonen map needs to be opti- 
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Fig. 1. — Schematic illustration of the structure (a) and func- 
tionality (b) of a Self-Organizing Map with / input neurons and 
M X N Kohonen neurons. The SOM visualizes the structure of the 
/—dimensional input space. In this case, the SOM illuminates a 
certain redshiftierror within the Kohonen map and as a function 
of the input space. 
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Fig. 2. — Accuracy (RMSE) versus regularization parameter 
value £ for the LRG-ELL data set (see Table 1). Different clas- 
sifications will result from different choices of the £ value. The 
regularization value is defined by the number of Kohonen neurons, 
which is maximum with respect to the training data set. The con- 
vex curve has a two local minima at £=4100 and £=5100. The 
roughness of this RMSE cost function shows that traditional gra- 
dient based optimization strategies, e.g. deterministic annealing, 
might result in sub-optimal solutions. Other methods, such as, 
genetic programming might find the global minimum much faster. 



mized for each data set. SOMs also allow non-experts to 
visualize the redshift estimates in relation to the multi- 
dimensional input space. This eliminates the often criti- 
cized "black box" problem of ANNs. As mentioned pre- 
viously, SOMs approximate the input feature space while 
ANNs typically separate them into sub-regions. Finally, 
SOMs are known to be powerful when very small data 
sets are available for training (see, Haykin 2009). 

4. CONCLUSION 
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Fig. 3. — Results from the three methods using SDSS u-g— r— i— 
z dcrcddcncd magnitudes as inputs for the SDSS DR7 Luminous 
Red Galaxies classified as ellipticals by the GalaxyZoo team. The 
bottom two plots show the SOM results for the two local minima 
described in Section 3 and shown in Figure 2 



SOMs offer a competitive choice in terms of low RMSE, 
algorithm comprehension (also see Goppert & Roscnsticl 
(1993)) and percentage of outliers. The final results are 
presented in Table 1 and plots for the LRG-ELL data 
set for the SOM, ANNz and GPR methods are shown in 
Figure 3 

As mentioned previously, obtaining the global mini- 
mum is important and, not surprisingly, can affect the re- 
sults. Figure 2 shows the two local minima (£=4100 and 
5100) listed for the LRG-ELL (Luminous Red Galaxies 
classified as Ellipticals by GalaxyZoo) data set in Table 
1. Clearly there are a number of other £- values and the 
RMSE will be greatly affected by the choice as seen on 
the y-axis of Figure 2 for a given £-value. Given these 
facts, the roughness of the RMSE cost function in Fig- 
ure 2 shows that traditional gradient based optimization 
strategies, e.g, deterministic annealing, might yield sub- 
optimal solutions. Other methods, such as, genetic pro- 
gramming might find the "global" minimum much faster, 



if a global minima exists with respect to the uncertainties 
of the RMSE. 

During completion of this manuscript another paper 
using SOMs for classification and photometric estima- 
tion was released (Geach 2011). Our work differs in that 
we mostly focus on a wider variety of low-redshift sam- 
ples drawn from the SDSS, while (Geach 2011) focuses 
more on the higher redshift samples akin to those used 
in Hildebrandt et al. (2010). 

We have shown that SOMs are a powerful tool for es- 
timating photometric redshifts and that with additional 
work they are sure to be useful in future surveys with 
limited numbers of follow-up spectroscopic redshifts. 
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